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Abstract. In traditional thermodynamics, temperature is a local quantity: a subsystem of 
a large thermal system is in a thermal state at the same temperature as the original system. 
Eor strongly interacting systems, however, the locality of temperature breaks down. We study 
the possibility of associating an effective thermal state to subsystems of infinite chains of 
interacting spin particles of arbitrary finite dimension. We study the effect of correlations and 
criticality in the definition of this effective thermal state and discuss the possible implications 
for the classical simulation of thermal quantum systems. 
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1. Introduction 

The question whether the standard notions of thermodynamics are still applicable in the 
quantum regime has experienced a renewed interest in the recent years. This refreshed 
motivation can be explained as the consequence of two successes. On the one hand, the 
spectacular progress of the experiments encompassed in the so called quantum simulators 
already allows for a direct observation of thermodynamic phenomena in many different 
quantum systems, such as ultra cold atoms in optical lattices, ion traps, superconductor 
qubits, etc [1, 2, 3, 4]. On the other hand, the inflow of ideas from quantum information 
theory provided significant insight into the thermodynamics of quantum systems [5, 6, 7, 8]. 
Specifically, qualitative improvements have been made in understanding how the methods 
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of statistical mechanics can be justified from quantum meehanies as its underlying theory 
[2, 6, 9, 10, 11, 12], 


One of the fundamental postulates of thermodynamics is the so called Zeroth Law: two bodies, 
each in thermodynamic equilibrium with a third system, are in equilibrium with eaeh other 
[13, 14], This is the law that stands behind the notion of temperature [13, 14]. In fact, 
the above formulation of the Zeroth law eonsists of three parts: (i) there exists a thermal 
equilibrium state whieh is eharaeterized by a single parameter ealled temperature, and isolated 
systems tend to this state [2, 6, 9, 10]; (ii) the temperature is local, namely, eaeh part of the 
whole is in a thermal state [13, 14]; and (iii) the temperature is an intensive quantity: if the 
whole is in equilibrium, all the parts have the same temperature [13, 14, 15, 16, 17, 18]. 
The last two points are usually derived from statistieal meehanies under the assumption of 
weakly interaeting systems. Nevertheless, when the interaetions present in the system are 
non-negligible, the points (ii) and (iii) need to be revised. Following the direction given 
by Refs. [15, 16], in this work we concentrate on the elarification and generalization of the 
aforementioned aspeets of the Zeroth Law of the thermodynamies for spin ehains with strong, 
short range interaetions. 

The general setting of the problem is as follows. The system with Hamiltonian H is in thermal 
equilibrium, deseribed by a eanonical state at inverse temperature (3, 





( 1 ) 


and we seek to understand the thermal properties of a finite part of the system. 

Obviously, in the presenee of strong interaetions, the redueed density matrix of a subsystem 
of the global system (espeeially in the quantum regime [18]) will not generally have the same 
form as (1). In lattiee systems, where the Hamiltonian is a sum of loeal terms interaeting 
aecording to some underlying graph, it is unelear how one ean locally assign temperature to 
a subsystem. More precisely, the redueed density matrix of the subsystem A (see Fig. 1) of a 
global thermal state is deseribed by 


PA = Tta{uj{H)) , 


( 2 ) 


whieh will not be thermal unless the partieles in A do not interaet with its environment A. 
Hence, given only a subsystem state pA and its Hamiltonian Ha, it is not possible to assign a 
temperature to it, sinee this would totally depend on the features of the environment and the 
interaetions that eouple the subsystem to it. 

In the eontext of quantum information, a first step to eireumvent the problem of assigning 
temperature to a subsystem was made in Ref. [15]. There, for harmonie lattiees it was shown 
that it is sufficient to extend the subsystem A by a boundary region B that, when traeed out, 
disregards the eorrelations and the boundary effects (see Fig. 1). If the size of such a boundary 
region is independent of the total system size, temperature ean still be said to be local. 

More explieitly, given a lattiee Hamiltonian H with a subsystem A, a shell region around it 
B and its environment C = {AVJ Bf, see Fig. 1, we aim to understand how the expeetation 
values of operators that act non-trivially only on A for the global thermal state uj{H) differ 
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Figure 1. Scheme of the subsystem of interest A, the boundary region B, and their 
environment C for a spin chain. Expectation values on A for the thermal state of the full 
Hamiltonian H (above) are expected to be approximated by expectation values for thermal 
state of the truncated Hamiltonian Hq (below) if the boundary region B is sufficiently large. 


from those taken for the thermal state of the truneated Hamiltonian Hab (with AB = AVJ B) 
as the width oi the boundary region B increases: 

1:t[Op'a]-Ti[Opa] < || 0 || oo /(®), ( 3 ) 

where pf = ArB^{HAB) is the state of A for the chain truncated to AB, and /(^s) is expected 
to be a monotonically decreasing function in £b- The width of the boundary region £b is 
defined as the graph-distance between the sets of vertices (regions) A and C. 

Surely, the differences (3) fully characterize the distance of pf from pA- Indeed, the trace 
distance, a distinguishability measure for quantum states, D{pf,pA) = \\\p'a ~ Pa\\i, has the 
following representation [19, 20]: 

D{p'a^ Pa) = ^max^Tr[0(p'^ - pa)] < /(®), (4) 

where / is the identity operator in the Hilbert space of A. 

In Ref. [17], it is proven that the correlations responsible for the distinguishability between 
the truncated and non-truncated thermal states are quantified by a generalized covariance. 
For any two operators O and O', full-rank quantum state p, and parameter r G [0,1], the 
generalized covariance is defined as 

covl{0,0') := Tr - Tr(pO) Tr(pO'), (5) 

and the average distinguishability of the two states measurements of some observable O can 
provide reads as 

A AO 

Ai[Ou{Ho)]-At[Ouj{H)] = 2 ds dt cov^//(i7,, O), (6) 

Jo Jo 

where Hj is the corresponding Hamiltonian term that couples B and C, Hq = H — Hj is the 
truncated Hamiltonian (see Fig. 1) and ojg = u{H{s) ) is the thermal state of the interpolating 
Hamiltonian H{s) := H — {1 — s) Hj. Hence, the generalized covariance is the quantity 
that measures the response in a local operator of perturbing a thermal state and ultimately at 
what length scales temperature can be defined. 

Temperature is known to be a local quantity in a high temperature regime. More specifically, 
in Ref. [17], it is shown that for any local Hamiltonian there is a threshold temperature 
(that only depends on the connectivity of the underlying graph) above which the generalized 








Locality of temperature in spin chains 


5 


covariance decays exponentially. Nevertheless, it is far from elear what oeeurs below the 
threshold, and, espeeially, at low temperatures (/3 1). Note that, in that ease, the right 

hand side of the truneation formula (6) eould be signifieantly different from zero sinee the 
integration runs up until 13/2, while the covarianee is expected to decay only algebraically for 
critical systems. 

In this work we show that, for one dimensional translation-invariant systems, temperature 
is loeal for any /?. Away from eritieality, we rigorously bound the truneation formula (6) 
by mapping the generalized eovarianee to the contraetion of a tensor network and exploiting 
some standard results in condensed matter. At eritieality, we use some results from conformal 
field theory [21, 22]. Finally, the results in [23], where the equivalenee of microcanonieal and 
eanonieal ensembles is proven for translation-invariant lattiees with short range interaetions, 
render our results valid also when, instead of being eanonieal, (1), the global state oj{H) 
is, e.g., mierocanonieal. The latter is defined as an equiprobable mixture of all the energy 
eigenstates in a narrow energy window (see [23] for details). 

In eondensed matter physies, this problem has been eonsidered in the eontext of 
approximating the expeetation values of infinite systems by finite ones, reeeiving the name 
of finite size scaling (see, e.g., [24] and referenees therein). Nevertheless, the finite-size- 
sealing methods are more focused to find the values for the critical exponents and the transition 
temperature by observing how measured quantities vary for different lattice sizes. 

The paper is organized as follows. In See. 2, the generalized eovarianee of ID system is 
mapped to the eontraetion of a 2D tensor network. In See. 3, we show that temperature 
is loeal at non-zero temperature {fi < oo) by identifying in the tensor network a gapped 
transfer matrix whieh leads to a elustering of correlations and ultimately to a elustering of 
the generalized eovarianee. In Sec. 4, locality of temperature is proven at zero temperature 
{(3 —)■ cxd) using different methods for gapped and gapless systems. While transfer matrix 
arguments work satisfaetory for gapped systems, eonformal field theory results have to be 
used at criticality. In Sec. 5, all our results are illustrated in detail for the Ising model, for 
which we study in addition the behaviour of the generalized covariance and compute explicitly 
the physieal distinguishability between the full and truncated Hamiltonians both at and off 
eritieality. Finally, we eonelude. 

2. Tensor network representation of the generalized covariance 

2.1. Mapping the partition function of a D-dimensional quantum model to the contraction of 
tensor network of D + 1 dimensions 

Let us eonsider a system of spins described by a short range Hamiltonian. The structure of 
the Hamiltonian is given by a graph G{V,E). The spins eorrespond to the set of vertiees V 
and the two-body interaetions to the edges E. Sueh a Hamiltonian ean be written as 



( 7 ) 


u£E 

where hu are the Hamiltonian terms aeting non-trivially on the adjaeent vertiees of u. 
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Figure 2. Diagrammatic representation of (a) the operator with H the Hamiltonian of a 
spin chain, (b) its decomposition (exp(—/3/miJ))™, (c) the tensor network that approximates 
performing the Trotter-Suzuki decomposition for a one dimensional short-ranged 
Hamiltonian and (d) the same tensor network after a convenient arrangement of the tensors. We 
use the Penrose notation; tensors are represented as geometric shapes, open legs represent their 
indices and legs connecting different tensors encode their contraction over the corresponding 
indices. 


In Ref. [26, 27], it is shown that, for any error e > 0, the matrix e of a loeal Hamiltonian 
can be approximated in one norm by its Trotter-Suzuki expansion. 


such that 



( 8 ) 




Ptn\\i < ^||e 


-I3H\ 


(9) 


where m > 360/3^ 177p/e and the products over u and v in Eq. (8) are realized in the same 
order. 

To illustrate the previous approximation, let us consider in detail the one dimensional case: 
a spin chain with nearest neighbour interactions. By decomposing in the standard way the 
Hamiltonian in its odd and even terms, the tensor network ptn becomes in this case 

( /y P P \ ^92- 

] , ( 10 ) 

where Hodd/even odd/even tind H TTodd T -^even- 

Let us think about each exp{/3/mhu) as a tensor. In this way, ptn can be seen as the 
contraction of several of such tensors, that is, a tensor network see Fig. 2 (c). Starting from a 
one dimensional quantum system, Ptn can be interpreted as a tensor network spanning two 
dimensions, with the extra dimension of length m. We will refer to this extra dimension as 
the (3 direction, while the original dimension will be called spatial direction. 

In Fig. 2 (a), the diagrammatic representation of Ptn is presented. Its tensors can be 
decomposed and arranged in order to form a square lattice of elementary tensors as shown 
in Fig. 2 (b). 
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Figure 3. Diagramatic representation of (a) the expectation value of a one site operator and (b) 
the generalized covariance (a two-point correlation function) between two one site operators. 
In both cases, the hnal result is computed as the ratio between the contraction of two tensor 
networks. 


2.2. Generalized covariance as the contraction of tensor networks 


The expectation value of a local operator is given by 

(0) = TiWi/)) = . (11) 

By using Eq. (9), the fact that ||e“^^||i = Z and some elementary algebra, the expectation 
value of a local operator A can be approximated by the ratio between the contraction of two 
tensor networks 


(O)- 


Tr(OpTN) 


<2||0||oo£. 


Tr(pT7v) 

This is represented diagramatically in Fig. 3(a). 
The generalized covariance can be rewritten as 

Tr f 

cov^4/s)(0,0') = 


Tr 


( 12 ) 


(13) 


where 0 = 0 — Ty{Ou{H)) for any operator O. Hence, in a similar way as it has been made 
for the expectation values, the generalized covariance can be also approximated as the ratio 
between two tensor network contractions as shown in Fig. 3(b) 


cov;;^4/3)(C’,o') 


^^{ORtnO Ptn 
Tt{ptn) 


<2||0||oo||0'||oo£. 


(14) 


From this perspective, the generalized covariance can be seen as a two point correlation 
function on a 2 dimensional lattice in which rm is the separation in the f3 direction and the 
distance between the non-trivial supports of O and O' is the separation in the spatial direction 
(see Fig. 3). 
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Figure 4. Diagramatic representation of the transfer matrix in the spatial direction (left) and 
the (3 direction (right). 


This construction can be generalized to approximate expectation values of local operators 
and n-point correlation functions of a D dimensional quantum model by the ratio of the 
contraction of two D + 1 dimensional tensor networks. 


2.3. Transfer matrices 


It is also very useful to define two extra objects: the transfer matrices along the spatial T and 
f3 directions T^. The first is obtained by contracting a column of the elementary tensors of the 
network, while the second is obtained by contracting several rows of elementary tensors, see 
Fig. 4. 

The number of rows that need to be contracted in order to obtain the transfer matrix in the 
(3 direction, Tp, is chosen such that its spectral gap between the largest and second-largest 
eigenvalues is independent of both f3 and m. This can be achieved by contracting m/(3 rows, 
leading to a transfer matrix with two largest eigenvalues pi and p 2 

(15) 




where A is the gap of the Hamiltonian. 


3. Locality of temperature at non-zero temperature 


Let us consider now the case in which (3 is of order one. The physical distinguishability in A 
between the full and the truncated Hamiltonians can be bounded by 

1 1 

Ti[Ou:{Ho)\—Ai[0 u:{H)]<l3 [ ds max cov^ (Hj,0)<2(3 [ dsmax max cov^ (hi, O), (16) 

J t6[0,1] " J i Te[ 0 ,l] 

0 0 

where Hj = hi + hR, i G {L, R}, and we have used the linearity of the generalized covariance 
with respect of its operators. 

Without loss of generality let us assume that the term of Hj that maximizes the generalized 
covariance is the one of the left, hi. Hence, the quantity to bound is cov^^(/ri, O). In order 
to do so, let us rewrite it as 


cov,. 


,, (iLlXsT^^YT^^TflR) iUlTsT^^YT^^TflR) { 1 l\X,T^^^+^T 41 r) 

fhL,U) = ------- 


(17) 
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Figure 5. Diagrammatic representation of a 3-point correlation function of the matrices Xg, 
Y and Tg for a system at zero temperature. The network is inhnite in both directions. 


where Zg := {1l\TsT‘^^‘^~^^Ts\1r) is the partition funetion, T is the transfer matrix in the 
spatial direetion, see Fig. 4 (left), and \Il/r) is the left/right dominant eigenveetor of T 
i. e. the eigenveetor assoeiated to its largest eigenvalue. The matrix Tg is the transfer matrix 
eorresponding to the boundaries BC where the elementary tensors of the network are different 
from the rest for s < 1 and Ti = T. The matrix Y eorresponds to the sliee of the region 
A where the operator O is supported, and the matrix Xg is the transfer matrix Tg with the 
insertion of the operators hi loeated at a distanee r/3 from O in the transverse direetion. 
The diagrammatie representations of the matriees Xg, Y and Tg are shown in Fig. 5. To 
simplify the ealeulations, the transfer matrix ean always be normalized sueh that its dominant 
eigenvalue is Ai = 1. 

To bound the generalized eovarianee (17) it is useful to rewrite it in terms of 2-point 
eorrelation funetions of the uniform system (s = 1) 


w. coM^B-,Xg,YT^-Tg) 

cov^fhL,0) = -- 

_ {Xg)T{Tg)T covTiis-, Tg, YT^^Tg) 

Zg Zg 

{Xg)T{YT^^Tg)T coYT{2iB + 1 ; Tg, Tg) 

Zg Zg 

_ {Tg)T{YT^^Tg)T covt ( 2 £b + 1 ; Tg) 

Zg Zg 

_ covTiis; Tg, YT^^Tg) covt ( 2 £b + 1 ; Tg) 

Zg Zg 

where 


(18) 


{X)t:={Il\X\1r), (19) 

covT{i-,X,Y) := {1l\XT^Y\1r) - {1l\X\1r){1l\Y\1r) , (20) 

and we have used that = {Tg)T and covt{^b] Xg,T^‘^~^^Tg) = covr(2£B -f 

1; Xg,Tg). 

In short range one dimensional systems, the absence of phase transitions at non-zero 
temperature [25] implies that the transfer matrix T is gapped, with a gap related to the spatial 
correlation length as 

e = -(ln|A2|)-'>0, 


( 21 ) 
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where A 2 is the second largest eigenvalue of the transfer matrix T. 

For gapped transfer matrices, the 2-point correlation function (20) can be proven to be upper- 
bounded by 

|covt(^;X,F)| < ||X|U)||||F|lij)||e-^/«. (22) 

The complete proof of the previous statement can be found in Lemma 1 of the Appendix A. 
Furthermore, Lemma 1 allows us to bound all the terms in Eq. (18), and, as it is shown in 
Appendix B, the following inequality holds for the left hand side of Eq. (16): 

Tr [Ouj{Hab)] - Tr [Ouj{H)] < 2/3c||/iiiUI|OiUe-^"/« + 0(e-2^^/«), (23) 

where c = 1 -f /Q^dscrL(s), with (Tl(s) = ((lL|T^Ti|lL) - \{Ts)T\‘^y^‘^ /\{Ts)t\ being the 
relative standard deviation of on the left dominant eigenstate |li). The quantity c is a 
constant of order one that depends on the model considered. Hence, the temperature is proven 
to be intensive for any one dimensional translationally invariant model at non-zero /). 


4. Absolute zero temperature 


4.1. Gapped systems 

Given a Hamiltonian with gap A, here we study the regime in which A. This implies 

that the lattice in its f3 direction is much larger than the correlation length (3 ^ ^ 13 , with 

-1 


I >n 


= A" 


(24) 


In the limit of temperature tending to zero, the 2D network that represents the partition 
function becomes infinite in the f3 direction (see Eig. 5). 

In order to see that the temperature is also local in this case, let us decompose the integral 
over t of the generalized covariance into two pieces 

/•l /•/3/2 

TvlOuiHAB)] - TriOuiH)] =2 ds dtcov^l^{Hi,0) 


JO JO 

pi pL pi pP/2 

= 2 ds dfcov^//(iTj,0)+ 2 / ds 

Jo Jo ^ Jo Jl 


dtcovl/f^{Hi,0) (25) 


where L is a cut-off that will be chosen afterwards to minimize a bound on the right hand side, 
and (3 will be made to tend to infinity. 

Concerning the integral over 0 < f < L, we will exploit the fact that the system is gapped, and 
hence its ground state is known to have a finite correlation length ^ in the spatial direction 
and to be represented by a Matrix Product State of bond dimenson D, with D oc poly(^) 
[40, 41, 42]. As argued in the previous section, a finite correlation length guarantees a gap in 
the transfer matrix in the corresponding direction. By performing an analogous calculation to 
the one described in the previous section, one obtains 

[ ds [ dfcov^/^(if7,0) < 2c||/i||oo^||0||ooe"^®/^ . 


( 26 ) 
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The second integral over t > L can be bounded by taking the transfer matrix in the (3 direction 
which is also gapped for gapped Hamiltonians. More specifically, the generalized covariance 
can be written as 

,.tl0(rr {GS\OTj,Hj\GS) {GS\OTl\GS) {GS\TlHfGS) 

^ ’ {GS\Tp\GS) {GS\Tj^\GS) {GS\Tj^\GS) ' ^ ’ 

where we have identified Tjs as the transfer matrix for which the ground state of the 
Hamiltonian \GS) is its dominant eigenvector. As previously, we make use of Lemma 1 
in the Appendix A and obtain 

< ||ff,||oo||Oi|ooe-*/«L (28) 

The integration is then bounded by 

hm / dt dscoy^Jf^{Hj,0)<^4HfU\0\U-G^^. (29) 

Jl Jo 

Putting the previous bounds together, and after an optimization over L, we get 

Tr [Ou{Hab)] - Tr [Ou;{H)] < 4\\0\U\Hi\\^ c C/j (^1 + y - In , (30) 

showing that temperature can be locally assigned to subsystems for arbitrarily large fJ and 
gapped Hamiltonians. 

4.2. Criticality 

A system at zero temperature is said to be critical when the gap between the energy ground 
state (space) and the first excited state closes to zero in the thermodynamic limit. The critical 
exponents zv control how the spectral gap A tends to zero 

AcxA-^% (31) 

where N is the system’s size and v is the critical exponent that controls the divergence of the 
correlation length 

e oc (32) 

The previous divergences are a signature of the scale invariance that the system experiences 
at criticality. If the critical exponent z = 1, there is a further symmetry enhancement and the 
system becomes conformal invariant. The group of conformal transformations includes, in 
addition to scale transformations, translations and rotations. 

In 1+1 dimensions, conformal symmetry completely dictates how correlation functions 
behave and how local expectation values of local observables of infinite systems differ from 
those taken for finite ones. Hence, conformal field theory establishes that 

Tt[Ou{Hab)] - Tt[Ou{H)] ~ ^ (33) 

up to higher order terms, where y is the scaling dimension of the operator Hj [38, 39]. If 
Hj is a standard Hamiltonian term, in the sense that the system is homogeneous, its leading 
scaling dimension is y = 2 . 

Once more, we see that by increasing the buffer region temperature can be arbitrarily well 
assigned. 
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Now we illustrate our results for the quantum Ising chain, which is described by the 
Hamiltonian 

1 h ^ 

= 2 - 2 

i=l i=l 

where a* and ex* correspond to the Pauli matrices, h characterizes the strength of the magnetic 
field and N is the number of spins. Notice that the interactions in the above Hamiltonian are of 
finite range, a crucial assumption in our derivations, see (7). This model has a quantum phase 
transition dXh = 1, so it well exemplifies the different regimes discussed above: criticality 
(only at zero temperature) and away from it (for zero and non-zero temperatures). 

5.1. Generalized Covariance 

First of all, as in the previous sections, we split the chain in three regions, which are shown 
in Fig. 6. For such a splitting, and in the context of Eq. (6), we compute the generalized 
covariance covL{f (O, O') taking for O a local operator in A, O = and for O' the 

boundary Hamiltonian between B and C, O' = Hj, given by 

= \ ® ■ (35) 



C BAB C 


Figure 6. Scheme of the subsystem A, the boundary region B and their environment C. The 
local operator acts on the subsystem A and the interaction term Hi corresponds to the 
red lines (connection between the subsystems AB and C). 


In order to compute cov^Ji^ {az^‘^iHi), we first diagonalize the Hamiltonian (34) using 
standard techniques from statistical mechanics, such as the Jordan-Wigner and the 
Bogoliubov transformation (see Appendix C). Once the Hamiltonian is diagonalized, we can 
straightforwardly construct the corresponding thermal state for every large but finite N, and 
compute covJ{f Hj) using expression (5). 

Figure 7 shows cov^Jf , Hj) as a function of t/(5 for several temperatures {j3 = 
5, 20,1000), and for h = 0.9,1 (i.e., near and at criticality). We take N = 40, which already 
describes well the thermodynamic limit (recall that we are only interested on the local state, 
and that the correlations decay exponentially). The area below the curves correspond to the 
first integral in (6), which measures how well the local state in A can be approximated by a 
thermal state in AB. 

The results in Fig. 7 are in agreement with properties (i) and (ii) from Lemma 2 in Appendix 
A. The first property implies that the covariance is symmetric with respect to f = /3/2, and it 
follows by taking I = t and n = jdin (A.8). Second, property (ii) implies that it is bounded by 
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a convex function of t with a maximum at f = 0 and t = (3 and with a minimum at t = (3/2. 
Therefore, the covariance satisfies the bound (A.9). 

On the other hand, the covariance is not monotonic in s (see Fig. 8). This is somehow 
counterintuitive, as it shows that the outcomes of two observables with no overlapping support 
(located in A and in the intersection between B and C) do not always become more correlated 
as s, which quantifies the strength of the interaction between B and C, increases. 


0.008 



0 0.5 to 0.5 1 

t/p tip 


Figure 7. Generalized covariance as a function of t for different values of s:s = 0,1/3, 2/3,1 
for the dotted, dashed, black and thick lines. The figures correspond to inverse temperature 
3 = 5 (top), P = 20 (at the middle) and /3 = 100 (bottom) and held strength h = 0.9 (left) 
and h = 1 (right). The grey area below the curves corresponds to the hrst integral of Eq. (6). 
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Figure 8. Generalized covariance as a function of s for different values of f: f//3 = 0,1 /3,1 /2 
for the thick, black and dashed lines. The figures correspond to inverse temperature /3 = 5 
(top), /3 = 20 (at the middle) and jS = 100 (bottom) and held strength h = 0.9 (left) and h = 1 
(right). The grey area below the curves corresponds to the second integral of Eq. (6). Notice 
that, due to the symmetry in t, the values f//3 = 2/3,1 are also considered. 


5.2. Locality of temperature in the quantum Ising chain 

In our analytical findings, the generalized covariance naturally appeared as a tool to solve 
the locality of temperature problem, see (6). This motivated the previous section, where we 
studied its properties in the context of the quantum Ising chain. Nevertheless, in order to 
obtain (6), one still needs to integrate covL{f ((9, if/) over s and r. While this approach is 
useful when dealing with arbitrary generic systems, here we are dealing with a specific model 
that is furthermore solvable, so we can take a more direct approach. Concretely, we first 


compute 





Pa — TrBc(i^(-f^)), 

with H = Hoc 

(36) 

and 





p'a = TrBc(i^(-f^AB)), 

with Hab = Hn, 

(37) 
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for different sizes N of the region AB. Seeondly, we measure the distinguishability between 
sueh states via the quantum fidelity, whieh is advantageous for eomputational reasons. For 
two states, pf and pA, the fidelity is defined as [19] 



(38) 


It satisfies 0 < F < 1 and F[pa,p'a] = 1 if and only if pA = pf- In order to relate this 
approaeh to our previous eonsiderations, we note the following relation between the traee 
distanee, D[pa, p^], and F(p^, pf), given in [20], 


l-F <D < s/I - F2. 


(39) 


Therefore, the fidelity provides us with upper and lower bounds to (4). In partieular, when 
D[pai p'^ —)■ 0 then F[p^, p'^] —)■ 1, and in that ease we say that the temperature is loeally 
well defined. 

From now on, we take for A a two spin subsystem, an infinite chain as the total system, and 
we compute F{pa, p'a) a function of the size of AB, with N = 2 + 21b, and the different 
parameters of the Hamiltonian. 

In order to compute pa and pf, it is convenient to apply the Jordan-Wigner transformation 
to (34), which maps spin operators crj. ^ ^ to fermionic operators a*, a| (see Appendix C for 
details). The Hamiltonian (34) takes then the form. 



(40) 


which is quadratic in terms of the fermionic operators. It follows that thermal states, as well 
as their local states, are gaussian operators. Therefore it is possible to describe them by their 
covariant matrix, whose size is only This allows us to compute pf in (37) for finite 

but large Ib', while in the limit N ^ oo, i.e. to compute pa in (36), we rely on the analytical 
results from [36]. The explicit calculations are done in Appendix C. 

5.2.1. Non-zero temperatures 

Figure 9 shows F{pa,p'a) as a function of fl and h, for N = 4 (left) and = 20 (right). 
Recall that N, with N = 2 + 21b, defines the size of the boundary region which is used to 
approximate pa by pf. Even if the boundary is small, = 4, the fidelity is close to 1 for all 
values of and h, and thus the temperature is locally well-defined. As expected, F(pa, p'a) 
increases with N (see Fig. 10). 

We also observe in Fig. 9 that the fidelity becomes minimal near h = 1, which is the phase 
transition point. As N increases, this minimum is shifted to h = 1. At this point the spatial 
correlations also increase, which suggests a relation between both quantities. 

In order to further explore this connection, we compute the scaling of F(pa, p'a) with N, and 
compare it to the decay of the correlations. The behaviour of F{pA, p'a) is plotted in Fig. 10, 
which clearly shows that the fidelity follows an exponential law with N, given by 



( 41 ) 
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where is a parameter that eharaeterizes the slope of the funetion. On the other hand, the 
eorrelations between a loeal observable in A, a*, and one in the interseetion of B and C, 



(42) 


can be obtained through the two-spin correlation function in [36]. Their asymptotic 

behaviour is also exponential with d. 



(43) 


where ^ is the correlation length. Now, identifying d, the distance between particles, with 
N/2, which is roughly the size of B, we obtain from the numerical results in Fig. 11 the 
following simple relation. 


e = 2^5. 


(44) 


Roughly speaking, the quality of the approximation pf is directly related to the strength of the 
correlations in the system. This relation is in good agreement with previous considerations in 
[43], where the correlation length is related to the error of the cluster approximation [43, 44]. 
In summary, temperature can be assigned to the local system for all h and non-zero fd by 
taking a small boundary region (with iV > 4, and thus Ib > 2). We have shown that this is 
directly connected to the exponential decay of the correlations with the distance, which makes 
the local state of a thermal state only be sensible to its closest boundary. 



Figure 9. Fidelity F{pA, p'f) as a function of /3 and the strength of the magnetic field h for 
N — 4 (left) and iV = 20 (right). The temperature is locally well-defined provided that 

F{pa,p'a) ~ 1 - 
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N 

Figure 10. Function of fidelity, — log(l — 
= 0.9 (right). The inverse temperature i 
dashed lines. 



F[pa, p'A\)y ^ function of iV for ft, = 1 (left) and 
s /3 = 5,20,100, 200 for the grey, black, thick and 


12 

9 

6 

3 

0 
0 

Figure 11. Correlation length, as a function of ft. The black spots correspond to the 
numerical values for 2 ^ 5 . The inverse temperature m (5 = h, 20, 75 for the black, dotted 
and dashed lines. 



5.2.2. Absolute zero temperature 

The same eonelusions apply at zero temperature, as the fidelity is also elose to 1 for all h and 
> 4. It also has a minimum near the eritieal point. 

Nonetheless, the sealing of the fidelity (or more preeisely 1 — F) with N ean differ from the 
previous ease. While the sealing is generally exponential at zero temperature, it beeomes a 
power law at the phase transition point (see Fig. 12), 

(45) 

This type of deeay is also obtained for the eorrelations as a funetion of the distanee, which 
again shows a direct connection between the quality of the approximation (quantified by 
F{pa-,p'a)) the strength of the correlations. 
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Figure 12. Function of fidelity, — log(l — F{pA, p'a))’ ^ function of N for /3 — oo. The 

field strength is ft. = 0.9,1. for the dashed and black lines. 


6. Conclusions 

In this work we studied the locality aspect of the Zeroth law of thermodynamics for quantum 
spin chains with strong but finite range interactions. Upon noting that in the presence of 
strong interactions the marginal states of a global thermal state do not take the canonical form 
themselves, we go on defining an effective thermal state for a subsystem. The latter being the 
reduced density matrix of the subsystem considered as a part of a slightly bigger, enveloping 
thermal system (see Fig. 1). Borrowing concepts from quantum information theory and 
employing methods from quantum statistical mechanics, we relate the accuracy with which 
the effective thermal state describes the actual state of the subsystem to the correlations 
present in the whole system (see Eqs. (4, 5, 6) and the discussion around them). We further 
utilize a Trotter approximation formula [26, 27] to build a tensor network representation of 
the corresponding states of the subsystem to provide upper bounds on the aforementioned 
accuracy, depending on the size of the enveloping thermal system, and such physical quantities 
as the spectral gap of the global hamiltonian and the temperature of the parent chain. At the 
quantum critical point, we use already existing asymptotical formulas from the conformal 
field theory. 

Lastly, we exemplify our analytical findings by analyzing the quantum Ising chain. The latter 
is complex enough to have a quantum phase transition point, but simple enough to allow 
for an exact diagonalization by standard tools of statistical mechanics, thereby serving as a 
perfect testbed for our analytical upper bounds. In particular, we find that, e.g., away from 
criticality, the envelope which is bigger than the system only by one layer of spins, is enough 
to approximate the actual state with a rather high precision (see, e.g.. Fig. 9). 

Our results for one dimensional systems with finite range interactions suggest that 
investigating the properties of the effective thermal states in higher dimensions and, possibly, 
harbouring long range of interactions, is an interesting direction for further research, which 
can have far-reaching implications in efficient simulation of the subsystems of large and 
strongly interacting quantum systems. Another interesting open question beyond the scope 
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of this work is whether these results can be generalized to other other types of equilibrium 
states, e. g. the so called Generalized Gibbs Ensemble and steady states of local Liouvillians. 
In a more practical vein, another field where our findings may find implications is quantum 
thermometry with non-negligible interactions [45]. 
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Appendix A. Proofs of the Lemmas 


In this section, we present the proofs of the lemma’s used in Secs. 3 and 4 to get statements 
on the locality of temperature for gapped systems. They consist of how different covariances 
decay for one dimensional systems with a gapped transfer matrix T. 

Lemma 1. [Infinite chain] Given a gapped transfer matrix T with eigenvalues \k labelled 
in decreasing order, i. e. |Afc| > \\ki\for all k < k', a right (left) dominant eigenvector |1^) 
((1l|A the largest eigenvalue Ai = 1, and a covariance between any two operators O and O' 
separated by a distance i defined as 

coy(i-0,0',T) = {1 l\O^T^O'\1r) - {1l\O^T^\1r){1l\T^O'\1r) . (A.l) 

Then, the covariance can be proven to decay exponentially in i 

|cov(f;0,0',T)| < \\0\U)\\\\0'\lR)\\e-^/^ (A.2) 

where ^ = — (In |A 2 |)~^ > 0 is the correlation length and the || |(^)|| = is the norm 

of the vector \(p). 


Proof Let us first introduce the new transfer matrix f = T — |lij)(l 2 ,| to rewrite the 
covariance as 


cov(£;0,0',T) = (U|Otf'0'|l^), 


(A.3) 


where we have used that — |l^)(l 2 ,|. By using the Cauchy-Schwarz inequality one 

gets 




1/2 


(A.4) 
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Let us now consider second factor separately. By inserting a resolution of the identity, a 
straight forward calculation leads to 

k >2 

< |A2p'5^{lL|0t|J;„){4H|0|li) < |A2p'||0|U)f (A.5) 

k >2 

where we have used that IA 2 I is an upper-bound for all the |Afc| with k > 2 and the Parseval 
inequality. 

Finally, we put everything together and get 

|cov(£;0,0',T)| < ||0|U)||||0'|l«)||e-'/« (A.6) 

where the 2nd largest eigenvalue | A 2 1 has been written in terms of the correlation length □ 


Lemma 2. [Periodic boundary conditions] Given a system with periodic boundary 
conditions, an Hermitian transfer matrix T with a gap A and a covariance between any two 
operators O and O' separated by a distance i defined as 


cov(f; n, O, O', T) = 


Tr (OT^O'T 




Tr(OT’^) Tr(0'T”) 


(A.7) 


Tr(T’") Tr(T’") Tr(T”) 

where 0 < f < n and n is the system size. Then, the covariance cov(f) = cov(f; n, O, O', T) 
as a function of i fulfills following properties: 

(i) Its real part is symmetric respect to the n/2 and the interchange of A and B, i. e. 

cov{n—i; n, O, O', T) = cov{i; n, O, O', T)* = cov(£; n, O, O', T) .(A.8) 

(ii) Given two operators O and O', there always exist two other operators Om and O'j^ such 
that 


|cov(£; n, O, O', T) | < cov(£; n, Om, O^, T) (A.9) 

and where cov(£; n, Om, O'j^^, T) is a convex function in i that is maximum at £ = 0 and 
1 , and reaches its minimum at £ = nl2. 


Proof. Statement (z) is a simple consequence of the following elementary equalities 

(Tr(OT^O'T’'-0)* = Tr ((OT^O'T'^-^O 

= Tr(T''-^0'T^0) = Tr(0'T^0T''-^). (A. 10) 

In order to prove {ii), let us focus on the first term in Eq. (A.7), since note that the second one 
does not depend on £. With this aim, we define 

Tr(Or^O'T"-0 
Tr(T") 


/(f;0,0')=3? 


(A.ll) 
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By introducing the transfer matrix in its speetral representation, f{i) ean be written as 




Tr(T”) 

^ ^ dkk' 

k,k'>2 


{iioiixi|o'|i)+y;c» 


k>2 


\k\ fXk 

IT niT 


t / \ \ n—l 


(A.12) 


where = 3? ((l|0|fc)(/c|0'|l)) and dky = 3^ {{k\0\k'){k'\0'\k)). Note now that 


Ai 


Ai 


Ai 


Ai 


n—l 


— + — = 2e cosh 


^-n |2 

ik 


where the eorrelation length ^kts defined as 

e—(I 


(A.13) 


(A.14) 


Note that as the eigenvalues of the transfer matrix are ordered, a larger k implies a shorter 
eorrelation length 

In a similar way, we can also simplify the terms in the last sum in Eq. (A.12). Note that 


Ai 


XiJ VAi 


k' 


n—l 




Ai.' 


XiJ VAi 


Ai 


_ jn _ n I i 

= e _j_ g Ik ^kk' 


= 2e cosh 


i-n/2 

^kk' 


.(A.15) 


where the length has been defined as • Puting the previous steps together, 

we get 




A? 


Tr(T’^ 


_ ti 

e Ck cosh 


k>2 


i-n/2 


^k 


E — ^ -—— 

dkk^ + 2 y, dkk'G cosh 


i-n/2 

^kk' 


. (A.16) 

k>2 2<k<k' A ^kk' J ^ 

Note that in general the covariance could oscillate in £, sinee and dkk' could take negative 
values for some k and k'. Nevertheless, given two operators A and B for whieh some 
Ck and dkk' are negative, there always exist two operators Om and sueh that their 
respeetive Cfc = \ck\ and dkk' = \dkk'\- For instanee, (^10^1^') = —{k\0\k') for the fc' 
and k-s with negative coefficents and {k\OM\k') = {k\OM\k') otherwise, and Ofj = O'. 
This eovariance f{i; Om, O'j^) is an upper bound to the absolute value of the previous one 
eovarianee f{i] 0,0'), 


f {i\ Om, Om) ~ 


A? 


E ti 

e ^^fe|ci;|cosh 


k>2 


Tr(T’ 

E _ Tl r -^ _ 

|4fc|e +2 2_^ \dkk'\e cosh 


i-n/2 


^k 


k>2 

>f{0,o,o'). 


2 <k<k' 


i-n/2 

^kk' 


(A. 17) 
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As the sum in Eq. (A. 17) is a linear eombination of eonvex funetions with positive eoeffieents, 
/(£; Omi O']^) is also eonvex. It is also obvious from the properties of the cosh() funetion, 
that f{i; Om, O'j^) reaehes its maximum at £ = 0 and n, and its minimum dXl = n/2. □ 


Appendix B. Some bounds on the generalized covariance 


In this appendix we bound the different terms in Eq. (18) to provide an upper-bound for the 
generalized covarianee used in the Section 3. 

Eet us start with the first term in Eq. (18). By using Eemma 1, we get 

|covr(£B;X„yT^«T,)| < ||X,|U)||||FT|T,|l«)||e-^"/«. (B.l) 

The coefficients can be bounded by 

||X,|U)||2 = (U|X,Xl|U) < ||/.,.||^(U|T,Tt|U) (B.2) 


||rT'-T,|l^)||2<||0||L||T'"+'T,|l«)||2 

< ||0||^ {\{Ts)t\^ + , (B.3) 

where we have used the structure of X^ and Y and the fact that Tr(Op) < ||0||ooTr(p) for 
any p > 0. It is also necessary to upper bound the inverse of the partition function. With that 
aim, let us write the partition function as 

Z, = {T,)| + (U|T, (T - Ts\1r) (BA) 

where we have used that (T — \Ir){Il\Y = T^-\ li?)(lL|. By proceeding similarly as in the 
proof of Eemma 1 , the partition function Zg can be lower bounded by 

|Z,| > \{Tg)T\^ - ||Tt|U)i|||T,|lR)||e-(2'"+')/«, (B.5) 


and its inverse is upper bounded by 


1-1 


<\{Tsh 


1-2 


1 + 


\{Ts)t\^ 


+ 




liTshr 

By putting the previous bounds together we get 

covTiiB;Xg,YT^^Tg) 




(B.6) 


< ||0||oo||/iL||oo(l+i^L(s))e-'«/«+0(e-2'«/«),(B.7) 

where (Tl(s) = {{1 l\TsT^\1l) — /\{Ts)t\ is the relative standard deviation of Tg 

on the left dominant eigenstate 11^) and we have omitted for simplicity the second order terms 
in 

In a similar way, the second term in Eq. (18) can be bounded by 

{Xg)T{Tg)T covriis; Tg, YT^^Tg) 


< ||/il||oo||0||oo(1 + cri(s))e 

+ C>(e“^^®/«), (B.8) 


where we have used that I{Xs)r/(^s)T| < ||^l| 
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The rest of terms in Eq. (18) ean be analogously bounded. Note that they will only eontribute 
to the seeond order. Putting everything together in Eq. (16), the physieal distinguishability on 
the region A between the truneated and untruneated thermal states is upper-bounded by 

Tr [Ou{Hab)] - Tr [Ou{H)] < 2Pc\\hL\U\0\\^e-^^/^ + 0(e-^^"/«), (B.9) 

where c = 1 + dscrL(s) is a eonstant of order one, depending on the model. 

Appendix C. Solving quantum Ising model 

In this appendix we find the states (36) and (37) using formalism of eovarianee matriees. 
Jordan-Wigner transformation 

Eet us first apply the Jordan-Wigner transformation, a* ® = (aj — ai)(aj+i + and 

(T* = ttiol — a\ai, to the Hamiltonian (34). We obtain, 

N ^ N 

Hn ^ ^ AijOiOj T — ^ ^ Oiof, (C.l) 

i,j=l ij=l 

with Aij = h 6 ij + and Bij = — (5ij+i) and where Oj and a\ denote 

annihilation and ereation operators, respectively. Erom this form of the Hamiltonian, we 
notice it is quadratic, and thus the thermal state (and their marginal states) are gaussian states. 
Therefore we can deal with them using the covariance matrix formalism. 


The correlation matrix 


In this formalism, we define the global correlation matrix, T, as 


r(x) = {xx^) = 


{OiOj) IIaTxTV 
(ajaj) IIa^xty 


II {aWj) IliVxTV 
II 


lAfxAf 


with X = 


ai 


On 

al 




, (C.2) 


where || ... \\nxn refers to a x iV matrix. Given T, we can obtain the correlation matrix 
corresponding to a reduced state by just selecting the corresponding matrix elements of T. Eor 
example, the correlation matrix of the fermions k,k + 1 is given by, 

^ (OfcCl^) (flkOik) (®fcOfc-|-l) 

(flk+ltltf) (Ufe+lU^_|_i) {pk+l(tk) (Ofc-i-iClfc+i) 

(44) (44+i) (4«fc) (4«fc+i) 

V (4-i-i4) (4-i-i4-i-i) (®I+i®fc+i) j 

Since the Jordan Wigner transformation is local, in the sense that it maps the k\h fermion to 
the fcth spin in the chain, this correlation matrix also corresponds to the two-spin subsystem 


/c,/c+l — 


(C.3) 
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at sites k and k + 1. This subsystem is preeisely the region of interest A in seetion 5.2, and 
thus (C.3) eorresponds to the eorrelation matrix of pa in (36). 

Given the redueed eorrelation matrix, the explicit form of pA can be easily obtained. As the 
reduced state of a thermal state is gaussian, there is a one to one connection between (C.3) 
and Pa- Indeed, for any gaussian state, with 




P = 


with M a coefficient matrix. 


Tr[e-x^MX^^ 

it is straightforward to prove that, provided that M is diagonalizable, 

r(A') = ‘ 


(1 + e-2^) 


(C.4) 


(C.5) 


or, equivalently, that 

M = -^iog(r(x)-i-i). 


(C.6) 


Explicit computation 

Now we explicitly compute (C.3) for a finite and an infinite chain, in order to obtain pf and 
Pa, respectively, using relation (C.6). 


• Finite chain 


For the case of a finite chain, we need to obtain the correlation matrix (C.IO) corre¬ 
sponding to the global state. It is then useful to first diagonalize the Hamiltonian (C.l) 
by applying the Bogoliubov transformation 


N 


bj ^ ^ T 'fjk)tli '4^jk)tf‘i, 


(C.7) 


k=l 


where f and f are real matrices and verify 4^% = J2k=i '4’% = 1- The Hamiltonian 
can then take the form, 

N 

H = Y,Ub{h-l/2), (C.8) 

k=l 


where are the fermionic excitation energies and bk and b\ denote annihilation and 
creation operators, respectively. The excitation energies, ^k, and the matrices f and f 
are obtained by solving the equation 


{A- B)cl) = 'fD, 


(C.9) 


where D is a diagonal matrix whose entries correspond to the excitation energies, ^k- 
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Once the Hamiltonian is diagonalized, it is easy to eompute the eorrelation matrix of a 
thermal state at inverse temperature (3 in the diagonalized basis, obtaining 


r(y) 


1 

l+e-fio 

OiVxTV 

OatxAT 

1 

1+ePO 


with Y 


where the non-zero matriees are diagonal. 




(C.IO) 


From that expression we can obtain the eorrelation matrix in the original basis, r(X), 
via 


r(X) = T^ViY)T, 


(C.ll) 


where T is the transformation matrix defined by the Bogoliubov transformation (C.7). 
That is, Y = TX, with 


and 


7 

p* 7* 


7 = ^(</> + V') and /i = -^(0 - '0). 


(C.12) 


(C.13) 


• Infinite ehain (iV —)■ cx)) 


In the ease of an infinite ehain, (C.3) ean be obtained relying on the analytieal results 
from [36]. The partial state of a two-spin subsystem is 


n^oo _ 

P2 — 




/ k- 






fc-l-i 


l=x^y^z 


, (C.14) 


where the average (cx^) and the two-spin eorrelation funetions {(crfcrf’'’^)}z={x,y, 2 } are 
given by [36]. In order to express the state in the fermionic basis, we ean eompute the 
redueed eorrelation matrix (C.3) from this state. 


with a 


/ 2(1 +a) -(/3 + 7) 0 

_ -(/5 + 7) 2(1 +a) /3-7 

0 /3-7 2(1-a) 

V -(/^- 7 ) 0 /3 + 7 

(a^), 13 = and 7 = (ajaj+i). 


-{(3-y) ^ 
0 

/3 + 7 

2(1 -a) y 


(C.15) 
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